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Spectral distribution functions of electron-phonon interaction a^F(Lj) obtained by ab initio 
linear-response calculations are used to describe various superconducting and transport properties 
in a number of elemental metals such as Al, Cu, Mo, Nb, Pb, Pd, Ta, and V. Their lattice dynamics 
and self-consistently screened electron-phonon coupling are evaluated within local density func- 
tional theory and using linear-muffin-tin-orbital basis set. We compare our theoretical a^F{Lj) 
with those deduced from the tunneling measurements and find a close agreement between them. 
Temperature dependent electrical and thermal resistivities as well as transport constants Xtr also 
agree well with the experimental data. The values of Xtr are close to the electron-phonon coupling 
parameter A. For the later a very good agrement with specific-heat measurements was found without 
any paramagnon contribution, except in Pd. We conclude that our method provides the description 
of electron-phonon interactions in tested materials with the accuracy 10%. 
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I. INTRODUCTION. 



Electron-phonon interaction (EPI) in metals is a sub- 
ject of intensive theoretical and experimental investiga- 
tions. The interest to this problem arises from a variety 
of physical phenomena such as electrical and thermal re- 
sistivities, renormalization of the electronic specific heat 
(electronic mass enhancement), and, of course, super- 
conductivity, for quantitative understandijHg of those a 
proper description of the EPI is requiredlll. Moreover, 
the reliable estimation of the EPI parameters in partic- 
ular case of high-temperature superconductivity may be 
decisive for recognizing the nature of this phenomenon. 
Unfortunately, even for some transition metals we have 
controversial experimental and theoretical data related 
to the estimation of the coupling constant A. Analy- 
sis is complicated by possible parallel processes of spin- 
fluctuations, for example, in the problemElof renormaliza- 
tion of specific heat and Tc or proximity effects in the tun- 
neling dataa. To extract the quantities under the interest 
one has to use some theoretical calculations and models. 
In this situation fully ab initio calculations of the one- 
electron spectra aiad phonon dispersions based on density 
functional thcoryQ (DFT) are most preferable. Applica- 
bility of the popular local density approximation^ (LDA) 
for the functional and treatment of the one-electron band 
structures as spectra of low-energy electronic excitations 
were checked many times and there, exist a reach theo- 
retical and experimental experienceO. This allows us to 
conclude that even having been formally ground-state 
theory, DFT is a good starting point for investigating 
the electron-phonon interaction. 



Many previous attempts to compute EPI, in partic- 
ular, for transition metals, have focused on calculat- 
ing merely the electronic contribution to this quantityQ, 
while the phonon frequencies ujq^ and the eigenvectors 
Vqi' were usually taken from inelastic neutron-scattering 
data. There, the self-consistent adjustment of the one- 
electron potential to the phonon distortion was replaced 
by either rigid-ion approximation (RIA)I3 or most popular 
rigid- mufhn-tin approximationcl (RMTA). For isotropic 
metals, having a large density of states at the Fermi en- 
ergy, the RMTA works well in many cases, since the ef- 
ficient electronic screening limits the change of the po- 
tential in the immediate vicinity of the displaced atom. 
However, there are known some problems of the RMTA 
in transition metals. For example, anisotropy of the mass 
enliajicement factor was not reproduced by the RMTA in 

Nm. 

Accurate ojq^ and rjqi, as well as self-consistently 
screened electron-phonon interaction can be calculated 
within the total^energy frozen-phonon approach using 
the supercellsEJ^EJ. However, there is a serious drawback 
of this method. Sufficiently large number of phonon wave 
vectors qmust be sampled in the Brillouin zone to get a 
good estimate of the average coupling strength A. A sep- 
arate frozen-phonon calculation is required for each q 
and for each studied mode. For small phonon wave vec- 
tors this requires very large supercells. With the crude 
sampling allowed by the limited size of the supercell, the 
accuracy of q integrated quantities like A is usually low. 

Another technique which can be employed for calculat- 
ing the self-consistfiiit change in the potential is the per- 
turbative approachllil applicable for any q . Key quantity 
of this method is an independent-particle polarizability 
function. After applying the first-order perturbation the- 
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ory and expanding the first-order changes in the one- 
electron wave functions over the basis of unperturbed 
Bloch states, the polarizabihty is expressed via the dotu 
ble sum over occupied and unoccupied states. Wintert£l 
has successfully applied this method to calculate A in 
Al. Unfortunately, the pcrturbative approach has sev- 
eral drawbacks. First, the slowly convergent sum over 
the excited states requires their preliminary calculation 
by diagonalizing the unperturbed hamiltonian matrix of 
very large dimension. Second, the self-consistency in this 
method is done by inverting the dielectric matrix of the 
crystal which is a relatively time-consuming problem. 

To date, the most efficient technique developed 
to calculate the lattice dynamics is the solid-state 
generalizationt^ of the Stcrnheimer methodO. This 
method is not limited to q's commensurate with the 
lattice as the frozen-phonon approach and it does not 
require the knowledge of all unperturbed states as the 
perturbative approach. The latter advantage can be 
achieved by constructing a rapidly convergent basis set 
for representing the first-order corrections to the wave 
functions other than the basis of unperturbed states. 
This is important because the first-order corrections as 
well as the unperturbed wave functions oscillate in the 
core region. In broad-band semiconductors and insula- 
tors this problem can be circumvented by using the pseu- 
dopotential method and most calculations of phonon dis.- 
persions performed so far use plane-wave basis setsE3. 
Unfortunately with decreasing band width the plane- 
wave expansion of the pseudo-wave-functions converges 
more slowly and it becomes less advantageous to use 
pseudopotentials. 

Recently an all-electron formulation of this—general- 
ized Stcrnheimer approach has been givenljE^. The 
first-order corrections are represented in terms of muffia-, 
tin (MT) basis set such as linear-mu0in-tin-orbitalsE3 
(LMTO) which greatly facilitates the treatment of lo- 
calized valence wave functions. For the first time, 
the method was shown to produce accurate phonon 
dispersions-jji—transition metals and transition-metal 
compoundsEjciJ. In this paper we present details of 
generalizing this method to compute the wave-vector- 
dependent electron-phonon couoling. (A brief report of 
this work has appeared alreadynj). We evaluate the spec- 
tral distribution functions a'^F{oj) of the EPI from the 
phonon linewi|dths according to the approach-devel- 
oped by AUencJ in the superconductivity theoryc3. The 
electron-phonon matrix elements are calculated in the 
LMTO representation. Due to incompleteness of the ba- 
sis sets in band-structure calculations the corrections to 
these matrix elements are shown to exist and explicitly 
taken into account. The incomplete-basis-set (IBS) cor- 
rections appear here in the same manner as in the calcu- 
lation p£the dynamical matrix within the linear-response 
theoryll3 or when calculating the forces within the total- 
energy £cpzen~phonon approach in terms of the LMTO 
methodN. 

We apply the developed scheme to compute electron- 



phonon coupling for a large number of elemental metals. 
We also present calculations of their phonon dispersion 
curves. The results of computed transport properties 
such as temperature-dependent phonon-limited electri- 
cal resistivities and thermal conductivities obtained as 
low-order variational solutions of the Boltzman equation 
are also given. The method of calculating the transport 
properties is analogous to-that used in the supeconduc- 
tivity theory and is basedcll on calculating the transport 
spectral functions a^j.F{uj). All the results presented in 
this paper are completely ab initio and no adjustable pa- 
rameters have been used in the calculations. 

The rest of the paper is organized as follows. In 
section II we derive the formulae for calculating the 
electron-phonon matrix elements and briefly review the 
method of finding superconducting and transport proper- 
ties. Section III presents the results of the calculations for 
phonon dispersions, electron-phonon interactions and re- 
lated properties for a number of elemental metals such as 
Al, Cu, Mo, Nb, Pb, Pd, Ta, and V. Section IV concludes 
the paper. 

II. METHOD. 

The central problem in calculating the electron-phonon 
interaction is the evaluation of changes in the electronic 
Hamiltonian caused by atomic displacements. This gen- 
erally requires the knowledge of the full low-energy exci- 
tation spectrum of the metal: the quasiparticle energies 
and the phonon frequencies. The calculations of vibra- 
tional properties are, in principle, within the scope of the 
density-functional based methods. Finding the quasipar- 
ticle excitation spectra is, on the other hand, much more 
difficult many-body problem. In the following we always 
assume that the quasiparticle energies are necessarily ap- 
proximated by the LDA energy bands. 

In the framework of the density-functional linear- 
response method the problem of calculating the phonon 
spectra and the electron-phonon interaction is reduced to 
finding the first-order variations in the one-electron wave 
functions, the charge density and the effective potential 
induced by the presence of a phonon with a given wave 
vector q. These quantities are connected by the so-called 
Sternheimer equation, which is the Schrodinger equation 
to linear order (we use atomic Rydberg units throughout 
the paper): 



(-V + Ve// - CkjJ^T^ h(- 



5^ 



5Rn 



SRn 



SRn 



)V'kj = 0, 

(1) 



where e^j , ipuj are the unperturbed energies and wave 
functions of vector k and band j; Veff is the effective po- 
tential of the DFT; d'^ekj/SR^,S'^iJkj/SRf„S^Veff/SRf, 
are the first-order changes in these quantities induced by 
the displacements of atoms in the positions R -I- 1 (R is 
the basis vector and t is the translation vector) moving 
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in the direction /i according to: Stu = SAiiexp{+iqt) + 
(5A^ea;p(— iqt) where SAn is the complex polarization 
vector. The sign "±" in the operator S^/SRf^ refers 
to the changes associated with the contributions to Stji 
which represent two travelling waves of vectors -l-q and 
-q. Note that (5±eki/<5i?^ = (kj|(5±V;///(5i?^|kj) van- 
ishes unless q=: 0. 

Equation must be solved self-consistently since 
the induced charge density S^p/SRfj, expressed via 
S^^pkj /SRfi, screens out the external perturbation 
S^Vext/SKfj,- The latter is simply given by the change in 
the bare Coulomb potential of the nuclei. We thus see 
that the scheme is the linearized analog of the original 
Kohn-Sham equations. If Eq.(|l|) is solved and the value 
of S^p/SRfj, is available, the phonon frequencies u^i, and 
the eigenvectors rjqi, (j/ numerates phonon branches) are 
obtained by diagonalizing the adiabatic dynamical ma- 
trix of the crystal. The latter is the second-order deriva- 
tive of the DFT total energy with respect to 6R^ and 
SR'^, and is conveniently represented by the Hellmann- 
Feynman theorem. In practice, however, there is a 
problem in using the Hellmann-Feynman expression con- 
nected with the fact that the unperturbed wave functions 
are obtained not as exact solutions of the Kohn-Sham 
equations but by expanding them into some basis set Xa ■ 



(2) 



Applying the Rayleigh-Ritz variational principle, this 
leads to finding the unknown coefficients from the 
matrix eigenvalue problem: 



v2 + K;/-6k,|xM^'-o 



(3) 



The consequence of the variational formulation is that, in 
general, the Hellmann-Feynman contribution is supple- 
mented with the correction expressed via the derivatives 
of the basis functions with respect to the displacements, 
(5*X^/(5i?^, which accounts for the incompleteness of the 
basis. The correction enters as the matrix element: 
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Ckjlkj) 



(4) 



and is obviously not equal to zero if the states f/'kj are 
only variationally accurate. The IBS corrections explic- 
itly disappear for the plane-wave basis set but must be 
taken into account for the basis of linear muffin-tin or- 
bitals which is used in our work. 

To treat properly the changes in the wave functions 
caused by variation of some external parameters one 
should check that the basis functions are well adjusted 
in the whole region of the parameters variations. Within 
the linear response, this is achieved by representing the 
first-order perturbations in ^/^kj in the form: 
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' SRu 



SRu 



(5) 



where S^A^i/SR,, are the changes in the expansion coef- 
ficients. This is advantageous comparing to the standard 
perturbation theory, where only the unperturbed states 
are used as a basis for representing S^ipi^j /SRfj, [This 
is described by only the first contribution in (||)]. The 
whole expansion (|5|) must be fastly convergent because 
the change S^x^/SR^^ in the LMTO basis is tailored to 
the perturbation just like the original LMTO basis x^ is 
tailored to the one-electron potential. 

In order to find the coefficients 6^ A^^ /SRfj, we first 
construct the functional of the dynamical matrix, min- 
imization of which with respect to 5^ipit.j /SRfj, leads 
to the equation (^, and then apply the Rayleigh-Ritz 
variational principle with the trial functions represented 
by Eq||. This leads to a set of matrix equations for 
S^A^^/6Rfj_. The detailed description of this method as 
well as its applications for calculating phonon dispersions 
can be found in Refs. [l9|-|2l|. 

In this paper we test the produced self-consistent 
change in the one-electron potential as the potential of 
electrons interacting with the phonon mode Wq^. Our 
task is to evaluate the electron-phonon matrix element 
ffk+qj' kj- T^^^ latter is conventionally written in the 
form: 



5k-|-qj',kj 



= {k + qj'\S'"'Veff\kj) 



(6) 



where both states tp^j and i/^k+qj' have the Fermi energy 
ep and where the change in the potential is transformed 
from the Cartesian system to the system associated with 
the eigenvectors rj^^ [RfJ') of a particular qiy - mode: 



eff 



1^ (AfflC.q,)l/2 ^ SR, 



eff 



(7) 



{Mji are the nuclei masses). 

It is not obvious but the expression for ^k+qj' kj 
should be corrected for the incompleteness of the basis 
functions. This fact immediately follows if we will inter- 
pret the matrix element (^) as a splitting of the degener- 
ate band ekj = Ck+qj' = e_F due to the phonon distortion. 
Consider first a non-degenerate state for which the first- 
order correction to the eigenvalue ekj is found as a change 
in the eigenvalue of the matrix problem (|^). It is given 
by 

5e^, ^ {kj\SVeff\kj) + 
(E4''^^^|-V2 + K//-ek,|kj) + 



{kj\-V' + V,ff-e^,\Y,Al^x^.) 



(8) 



where 5x^ and SVeff denote the changes in the basis 
functions and the potential due to some change in the 
external parameters of the hamiltonian. This formula 
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contains both the expression of the standard perturba- 
tion theory (first term here) as well as incomplete-basis- 
set corrections [second and third contributions in (||)]. 
It would be advantageous to use this formula because 
the eigenvalues of the matrix problem are variation- 
ally accurate for the whole range of parameters variation. 
The same valid if we develop a perturbation theory for 
the change in the degenerate band e^j — fk+qj' within 
(^. This leads to the result: 



(9) 



where denotes the variation of the MT basis func- 

tions due to the phonon distortion of the qi' -mode and is 
connected with the variation S'^^^jSR^ in the same way 
as for the induced potential, Eq. (|^). The last two con- 
tributions in represent the IBS-corrections which are 
not vanished unless i/jki, V'k-i-cn' are the exact solutions. 

Formally, the expression (0) can also be obtained by 
repeating the standard quantum-mechanical derivation 
of the Fermi " golden rule" for the wave functions repre- 
sented in terms of the basis set according to (||). The 
matrix element of the electron-phonon interaction is ob- 
tained by considering the scattering rate for transitions 
from an initial, unperturbed state V'kj(i) into a final per- 
turbed state i/'k'j' {t) at the time moment t which is given 
by the overlap integral squared: 



l(kjW|k'/W)|^ 



(10) 



Since the final state corresponds to the displaced lattice, 
the best variational estimate for i/'k'j' (0 must include 
the orbitals centered at the new atomic positions and ad- 
justed to the new one-electron potential. To linear-order 
with respect to the displacements this leads to finding 
the change in the basis Xq which gives rise to the IBS- 
corrections entered (^. 

The expression (^ is thus the linear-response analogy 
of evaluating the clcctron-phonon matrix elements via 
the splitting of the bands in the frozen-phonon super- 
cell method as done in Ref. ^ It is less sensitive to 
the errors in the wave functions introduced by the varia- 
tional principle, has a correct long-wavelength behavior 
and allows one to avoid the inclusion of rf — / transitions 
in d— electron systems. 

For the electron-phonon spectral distribution func- 
tions a'^F{ijj) we employ the expressiono in terms of the 
phonon linewidths 7qj/ 



E— 



27riV(eF) ^ uj. 



(11) 



where N{eF) is the electronic density of states per atom 
and per spin at the Fermi level. When the energy bands 



around the Fermi level are linear in the range of phonon 
energies, the linewidth is given by the Fermi "golden 
rule" and is written as follows 

7q^ = 27rwq^ E l5k+ci,',kjP^(<^kj - ei.)(5(ek+qj' - ep)- 

(12) 



The spectral distribution function (11) and its first recip- 



rocal moment A are usually used to describe such impor- 
tant manifestation of the EPI as superconductivity and 
some normal-state properties. One of such properties is 
an enhancement of the electronic mass for the electron 
at the Fermi energy when its velocity Wk is reduced by 
the factor 1 -f Ak due to the interaction with phonons. 
The value of Ak is given by the reciprocal moment of the 
so-called k— dependent electron-phonon spectral function 
a'^F{uj). This renormalization is observed in the de Haas- 
van Alphen and cyclotron-resonance experiments. As its 
consequence, the low-temperature electronic specific heat 
is also renormalized. For the latter effect it is sufficient 
to know only the Fermi-surface averaged value of Ak, i.e. 
A. 

Full description of the superconducting state can be 
obtained by solving the Eliashberg gap equations which 
relate the energy-gap function and the renormalization 
parameter for superconducting state to the electron- 
phonon aud the electron-electron interactions in the nor- 
mal stateE3. The electron-phonon-coupling function is 
given by a^F{uj). The Coulomb interaction is usually 
represented by some constant /i* . Detailed nature of the 
effective Coulomb repulsion is not very well known. For- 
tunately, various definitions for fi* have only a weak in- 
fluence on the solution of the gap equations and its values 
can, e.g., be found by adjusting the calculated transition 
temperatures to the experimental ones. 

Electron-phonon scattering has a dominant contribu- 
tion to the electrical resistivity for reasonably pure met- 
als except very low-temperature region where the im- 
purity and electron-electron scattering are important. 
The influence of the EPI on the transport properties ape 
described in terms of the transport spectral functiontJ 



Ficu) 



at 



1 NT^ 



ai^F{uj) where: 



l5k;-1^-|V(k)t;.(k('))x 



NUf){vI) ^ ^ 
^(eki - eF)'5(ek'i' " eF)^{^ - t^k'-ki.) 



(13) 



Here (u^) is the average square of the x component of the 
Fermi velocity. In the lowest-order variational approxi- 
mation (LOVA) for the solution of the Boltzman equation 
the expressions for electrical and thermal resistivities are 



N{ep){v'^) J uj sinK^x 





al.F{Lo), 



(14) 
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w{T) 



i:kBN{ep){v'^) J uj sinh'^x 




duj x'^ 



Ax 



2x2 



[aiF{uj) + —at,,tF{oj) + — a^„F(^)] (15) 



with X = uj/2kBT. The LOVA resuhs (jlj), (|lj) give the 
upper bound to the resistivities and aUow us to test the 
calculated a'fj.F{iLj). 

III. RESULTS. 



a. Technicalities. 

Our calculations of phonon dispersions and electron- 
phonon interactions for the elemental metals such as fcc- 
Al, Cu, Pb, Pd, and bcc- Mo, Nb, Ta, V are performecUii 
the framework of the linear-response LMTO methodcil. 
The details of the calculations are the following: We 
find the dynamical matrix and the phonon linewidths for 
these materials as a function of wave vector for a set of ir- 
reducible q-points at the (8, 8, 8)-reciprocal lattice grid 
[29 points per l/48th part of the Brillouin zone (B/ 
The (/, J, K) reciprocal lattice (or Monkhorst-Pack 



grid is defined in a usual manner: q^fc = jGi 

k 



-^Ga, where Gi,G2,G3 are the primitive translations 
in the reciprocal space. The self-consistent calculations 
performed for every wave vector involve the following 
parameters: We use 3k — spd— LMTO basis set (27 or- 
bitals) with the one-center expansions performed inside 
the MT-spheres up to Imax = 6. In the interstitial region 
the basis functions are expanded in plane waves up to the 
cutoff approximately corresponding to 70, 140, and 200 
plane waves per s,p, and d— orbitals, respectively. All 
semicore states lying higher than —ARy are treated as 
valence states in separate energy windows. The induced 
charge densities and the screened potentials are repre- 
sented inside the MT-spheres by spherical harmonics up 
to Imax = 6 and in the interstitial region by plane waves 
with the cutoff corresponding to the (16, 16, 16) fast- 
Fourier-transform grid in the unit-cell of direct space. 

The k-space integration needed for constructing the 
induced charge density and the dynamical matrix is per- 
formed over the (16, 16, 16) grid (145 points per l/48th 
part of the BZ) which is twice denser than the grid of the 
phonon wave vectors q. We use the improved tetrahedron 
method of Ref. However, the integration weights for 
the k— points at this (16, 16, 16) grid have been found 
to take precisely into account the effects arising from 
the Fermi surface and the energy bands. This is done 
with help of the energies ekj generated by the original 
full-potential LMTO method at the (32, 32, 32) grid (897 
points per 1/48 BZ). The procedure is in details explained 
in Ref. |l] and allows us to obtain more convergent results 
with respect to the number of k points. 

The k-space integration for the phonon linewidths 
7qi/ is very slowly convergent because it involves two 



5 functions according to Eq. (|12|). It is performed 
with help of the (32, 32, 32) grid in the BZ by means 
of the tetrahedron method of Ref. ^ The largest 
numerical error of a'^F{uj) comes from the integration 
over q in the expression (pT|). Its magnitude, we es- 
timated by performing the integration over merely the 
band-structure factor [which is 7qi/ approximated by 
Skij' <5(ekj-e_F)'5(ek+qj'-eF)] using respectively (8, 8, 8) 
and (32, 32, 32) grids and found to be not larger than 7% 
in all cases. 

A few words should be said about the lattice parame- 
ters used in the calculations. It is known that the equi- 
librium cell volume V found theoretically by the cor- 
responding LDA-based total-energy calculation is fre- 
quently obtained slightly lower than the experimental 
volume Vq. This usually leads to the calculated at Vb 
phonon frequencies which are softer comparing to the 
experimental ones. Often, a better agreement with the 
experiment can be obtained by performing the linear- 
response calculations at the theoretical volume. This is, 
in principle, a justified procedure from theoretical point 
of view. Unfortunately, the prescription does not work 
when calculating the phonon linewidths and af^F(w) : 
the results of calculated electrical and thermal resistiv- 
ities agree less well with the experiment. The reason 
for this discrepancy is connected with the sensitivity of 
these quantities to the shape of the Fermi surface. It 
turns out that the use of the Fermi surfaces calculated 
at the experimental lattice constants considerably im- 
proves the results. We can thus use theoretical volumes 
in the linear-response calculations of phonon dispersions 
and the electron-phonon matrix elements. To find the 
phonon linewidths and a^^F^uj) we can use, on the other 
hand, the energy bands entered ( p^ ) which are generated 
at the experimental lattice constants. We understand 
that it is not well justified procedure to use different lat- 
tice parameters in one calculation, but it somewhat helps 
to minimize the errors connected with the LDA by simple 
means. The actual volume ratios V/Vq used in our calcu- 
lations to find the changes in the one-electron potentials 
are listed in Table I. 

Another comment concerns the choice of the exchange- 
correlation potential. The general strategy employed by 
us is to use the exchange-correlation formula which gives 
the best prediction of the cell volumes. The Barth- 
Hedin-like formula after Ref. ^ is employed for all the 
metals except Cu and V. For the 3(i-metals we have 
found that this formula gives the theoretical volumes 
which are too small ^/Vq ~ 0.9). For Cu and V the 
Ceperley-Alder forn£3 of the exchange-correlation po- 
tential parametrized after Ref. 33 is used which gives the 
ratios V/Vq being much closer to unity (see Table J). 

Finally note that in the previous publicationsEauJa 
we have used a different method for treating the full- 
potential terms in the calculation which was based on 
the atomicj-eells and the one-center spherical-harmonic 
expansion's. This method is not directly applicable 
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to calculate phonon dispersions for materials with open 
structures such, e.g., the diamond structure and requires 
the replacement of empty sites of the lattice by empty 
spheres. This complicates the evaljiation of the dynam- 
ical matrix. In recent publicationEll we have employed 
another approach based on the plane-wave expansions 
for the LMTOs in the interstitial region and have ap- 
plied the method to calculate the phonon spectra in Si 
and NbC. While the materials considered in this work 
have close-packed bcc or fee structures we also apply 
this new method which is more general for practical 
use. Somc-n£,thj3 results for Al, Nb, and Mo previously 
publishedtj'c3nJ do not noticeably differ from those pre- 
sented below in this paper. 

b. Lattice-dynamical properties. 

In Table I we report the values of the calculated phonon 
frequencies at the high-symmetry points X and L for the 
fee metals Al, Cu, Pb and Pd as well as at the points H 
and N for the bcc metals AIo, Nb, Ta, and V . For the 
comparison, the experimental frequenciesEil are also listed 
in Table I along with the theoretical-to-experimental vol- 
ume ratios which have been used in the calculations. 

Our results for the phonon dispersions along several 
symmetry directions together with the corresponding 
densities of states for these materials are displayed in 
Fig.l(a)-(h). The theoretical lines result from the in- 
terpolation between the calculated frequencies which are 
denoted by circles. r— Many neutron-diffraction measure- 
ments are availableEil for nearly all the metals considered 
here and these data are also shown in Fig.l by trian- 
gles. The only exception is V for which the dispersion 
relations cannot be studied with neutrons since V is an 
almost totally incoherent neutron scatterer. While some. 
X-ray diffraction measurements exist in the literatureEfl 
their accuracy seems to be less satisfactory than corre- 
sponding neutron-scattering data for other materials and 
we do not show the experimental points for V. 

From Fig.l we see that the agreement between theory 
and experiment is good. Most of the calculated frequen- 
cies agree within a few percent with those measured. This 
also follows from the numerical values listed in Table I. 
In particular, for Al [Fig. 1(a)] a very good agreement 
is found in all the directions. As we have mentioned 
already, this calculation is performed at the theoretical 
volume {V/Vq = 0.955). We have also checked the set-up 
with the experimental volume and found a considerable 
softening (about 20%) of the transverse modes. This il- 
lustrates the importance of performing lattice-dynamical 
calculations at the theoretical volumes. 

The most important consequence of our calculation for 
lead [Fig. 1(b)] is that the pronounced dip of both the 
longitudinal and transverse branches near the X-point 
is well reproduced . We have also found a slight over- 
estimate of the transverse phonon frequencies near this 
zone boundary which can be attributed neither to the dis- 



crepancy in the cell volume {V/Vq = 1.002) nor to the ne- 
glection of the semicore states. (We have in fact included 
both 5d, and 6d-states in the main valence panel). This 
kind of disagreement has also been recently reported in 
Ref. ^ using the linear-response pseudopotential tech- 
nique. It is possibly connected with the use of the local 
density approximation or the lack of spin-orbit coupling 
effects in our calculation. 

The most interesting cases are V, Nb, and Ta 
[Fig.l(c),(d),(e), respectively]. The materials belong to 
Group V of the periodic table and all of them show 
anomalious behavior of the phonon dispersion curves. 
The presence of anomalies is, first of all, connected with 
the well-known dip of the longitudinal mode in the (00^) 
direction. The dip is correctly reproduced by our cal- 
culation. Another important features of our calculation 
are (i) the softening of the transverse mode along the 
(00^) direction at long wavelengths which is rather sharp 
in both V, Nb and is weaker in Ta as well as (ii) the 
crossover of two transverse branches in the (f ^0) direc- 
tion. The latter is in fact predicted for Ta because the 
measured dispersions in this direction are absent for the 
Ti-branch except the zone-boundary point N. 

The theoretical phonon dispersions for Mo also agree 
well with the experiment. The consequences of our cal- 
culations here are the reproduced softening near the H 
point and the absence of the large dip along the {£,£,£,) 
direction near ^ ~ 0.7. The dip is presented in the dis- 
persion curves of nearly all bcc metals except of Mo 
and Cr, and is certainly the feature of the behavior of 
the crystalline structure factor which enters the dynam- 
ical matrix as a sum over lattice vectors with the phase 
shift exp{iqt). Its absence indicates a considerable wave- 
vector dependence of the electron-phonon matrix ele- 
ments which seems to be well reproduced by our method. 

Finally we compare the results of our calculations for 
Cu and Pd which are presented in Fig. 1(f) and Fig. 1(g). 
The dispersion relations for these materials do not show 
the anomalies and are smooth. A slight overestimation 
of the theoretical phonon frequencies is found for both 
these metals . The overestimation can likely be corrected 
by the use of the cell volumes in the calculation which 
would be slightly larger than those found within the 
LDA. However, such a procedure is not justified theoret- 
ically. We think tiiat employing the gradient-corrected 
density functional known to predict much better the 
equilibrium lattice parameter will allow us to improve 
these results. 

c. Superconducting properties. 

We now discuss applications of our ab initio linear- 
response method to calculate the electron-phonon cou- 
pling and superconducting-state properties. First, we 
present our results for the spectral distribution functions 
and A. Second, we describe our applications to solving 
the Eliashberg gap equations. 
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Calculated a'^F{uj) for Al is shown in Fig. 1(a) (full 
line). The positions of the maxima here are conditioned 
by the form of the phonon density of states with the 
low-frequency phonon peak suppressed by the coupling 
function a^{uj) (dashed line). The broad phonon spec- 
trum in Al is extended up to the maximal frequency 
i^max ~ 470i^. The theoretical a'^F{uj) in Fig. 1(a) is corrtj 
pared with the results of the tunneling measurementsu 
(squares). We find a rather good agreement between the 
two curves. In fact, our a'^F{uj) is also found to be prac- 
tically identical to the empirical pseudopotential result 
of Ref. ^ based on the rigid-ion approximation. The 
latter is known to work well in simple metals. General 
agreement is found between our and|-tlie ab initio frozen- 
phonon results of Dacorogna et. altn for the dispersion 
of the phonon linewidths along the high-symmetry di- 
rections. The only exception is that, in the (00^) direc- 
tion, our longitudinal branch of 7qi, exceeds theirs by a 
factor of 2. This is presumably connected with replac- 

fthe S functions in (|l2|) by Gaussians used in Ref. 
However, the relative weight of our high 7 values 
in the integrated quantities, such as a'^F{uj) and A, is 
found to be very small. Our value of A is 0.44 which 
is very close to the value of pXt„„ = 0.42 extracted-J'rom 
the tunnelling measurementsB. The frozeii-phononll3 and 
linear-response calculations of WinterEj gave, respec- 
tively, A — 0.45 and 0.38. The value of Xs~h extracted 
from the electronic specific-heat coefficient 7 and our cal- 
culated density of states N{eF) using the relation 



1 + A,-,, = 



37 



(16) 



is 0.43 (see,|-Table II). In order to check previous 
conclusionsESS about the inapplicability of the RMTA 
for sp metals, we also performed such a calculation and 
indeed found Xrmta = 0.14. 

Lead is a well-studied classical example of strong- 
coupled superconductor with =^ 7.19K and its tun- 
nelling spectra have been studied a long time agdS. Ob- 
tained a'^F{uj) using our linear-response method is pre- 
sented in Fig. 2(b) where it is compared with the results 
of the measurementsEa. The two curves are similar. Our 
calculated A — 1.68 is found to be 8% larger than the tun- 
nelling value 1.55 and only 2% larger than the value 1.64 
extracted from specific-heat data (see Table II). This dis- 
agreement is well within the accuracy of our calculation. 

We now report our results for V, Nb and Ta which 
are the best-studied elemental superconductors because 
of their relatively high-Tc values. Especially, for Nb 
which has the highest Tc ~ 9.25K among the elemen- 
tal metals, there exist maajyi-epcperimental investigations 
of the tunneling SDectracJES and theoretical RMTA- 
based calculationsc3if3. Unfortunately, some of the re- 
sults which have been reported in the literature are con- 
troversial. First, the RMTA calculations give the values 
of A varying from 1.12 to 1.86. Second, the tunneling 
estimates of the coupling constant for Nb and V are dif- 



ficult because of the oxidation of the surface layers. Hav- 
ing lower transition temperature, such oxides act on the 
tunneling spectra due to the proximity effect. The efc 
perimental A in Nb varied in the past from the valuesLj 
0.58 with negative or anomaliously small fi* to the 

valudlj 0.9. At present, a satisfactory explanation of the 
anomalious behavior of the thermally oxidized tunnelling 
junctions in Nb appears to be possible which gives the 
valueMy of Xtun = 0.92 ~ 1.22. 

In Fig. 2(c), (d),(e) we present our calculated a'^F{uj) 
(full lines) for V, Nb and Ta, respectively. They all 
have rather broad spectra extended up to ujmax ~ 
370^,3107^, and 240i^. The coupling function a^{ijj) 
(dashed line) in these metals only slightly deviates from 
constant in major part of the frequencies. The approx- 
imation — const, works very well in Nb and Ta. 
This qualitative result implies that the electron-phonon 
coupling can be factorized into electronic and phonon- 
dependent factoralJO. 

In Fig. 2 the calculated spectral functions are compared 
with the results of the tunneling measurements (squares) . 
As it can be seen for l/|-[Fig.2(c)], our a^F(w) disagrees 
with the measured oneca because of the appearance of 
the upper phonon peak not presented in the experiment. 
Even though the theoretical a'^F{u)) should be broad- 
ened because the 5 function in Eq. (jlll) ought to be a 
Lorentzian of half- width 7qi,, the electron-phonon cou- 
pling estimated by us (A = 1.19) is 40—50% stronger than 
the obtained Xtwn — 0.82. The same situation is found for 
Nb. Our calculation here [Fig. 2(d)] also does not show the 
suppression of the longitudinal peak. The latter is absent 
in nearly all the exjjeriments for this metalcil. [A typical 
measured spectrumQ is shown in Fig. 2(d) by squares]. 
As a result, |-the calculated A — 1.26 is 20% higher than 
Xtun = 1.04Q. The discrepancy found by us has alijeady 
been reported in the past RMTA-based calculationsC3E3. 
To check the consistency of our results with the earlier 
ones, we have performed our own RMTA calculations 
and obtained complete agreement between them. We 
thus conclude that the full inclusion of screening does 
not resolve the problem of the suppressed longitudinal 
peak. 

Unfortunately, our comparison with the experiment is 
complicated by the proximity effect and the extraction 
of the tunneling densities of states depends on the way 
how the measured data are processed. For example, in 
Nb the value of Xtun = 1-22 deduced from the tunneling 
experiments (which is only 3% lower thari that found by 
us) has been reported in the literaturec3. The obtained 
a^Fiuj) [denoted in Fig. 2(d) by triangles] is found much 
closer to our calculation. 

A better understanding of the present situation can 
be achieved_tw comparing the theoretical and the 
experimentaSB tunneling spectra for Ta since its su- 
perconducting properties are close to those of V and Nb 
but this metal is much less reactive with oxygen. Such 
a comparison is given in Fig. 2(e). We find rather good 
agreement between the both curves. In particular, the 
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upper phonon peak is not suppressed in the measured 
a'^F{uj) and its amplitude is comparable with that calcu- 
lated by us. As a result, the theoretical A = 0.86 agrees 
within 10% with the At„„ = 0.78. 

From view of the data on Ta we cannot consider the 
discrepancy found for V and Nb as a drawback of ei- 
ther our linear-response method or the use of the local 
density approximation. Partially this conclusion is also 
verified by alternative estimates of the coupling constant 
based on the specific-heat data and the de Haas-van 
Alphen (dHvA) experiments, but it should be noted that 
both cyclotron masses and the specific-heat coefficient 
are also enhanced by the electron-electron interactions. 
Evaluation of the average coupling from the specific-heat 
measurement£j yields an enhancement of 1.00 for V and 
1.17 £qx Nb (see Table II). If one uses a more recent 
valueEJ of 7 for V rather than listed in Ref. ^ one obtains 
the enhancement Xs-h = 1-17 which is nearly coincides 
in this metal with our A ~ 1.19 . Comparison between 
the I|JQA band masses with those measured by dHvA 
effectEil yields the enhancement of 1.33 for Nb, which is 
close to the value 1.26 found in our calculation. Another 
important result is that the measured variation of the 
mass enhancement for the various cyclotron orbits in Nb 
agrees also well with our calculation. Namely, we have 
found a decomposition of A by the Fermi-surface sheets: 
(i) octahedron, (ii) jungle gym and (iii) ellipsoids, and 
have obtained the contributions A' = 1.44, A" = 1.37, 
and A"' = 1.08. They can be compared with the mea- 



sured A 



dHvA 



1-71,A:J^„^ = 1.43, and A 



dHvA 



1.10. 

(Note that the calculated within the RMTA anisotropy 
of tbe, mass enhancement strongly disagrees with these 
dataEd.) Moreover, we have estimated the transport con- 
stants Xtr for these metals both from the calculated and 
measured resistivity data. The values of Xtr are usually 
believed to be close to the superconducting ones. (A com- 
plete report of our calculated transport properties will be 
given in the following subsection) . For V, the value of Af ,■ 
found by us is 1.15 and, for Nb, Xtr — 1.17. The calcu- 
lated electrical and thermal resistivities are also close to 
those measured. It therefore seems that our electron- 
phonon coupling is accurate (within the computational 
accuracy of order 10%) while the effect of the electron- 
electron interactions is small in these metals. 

As the next two examples, we report the results of 
our applications for Mo and Cu. There are no tunnel- 
ing data for these materials because of their low Tc and 
the weakness of phonon effects. The calculated spectral 
functions are presented in Fig. 2(f), (g). Both curves qual- 
itatively agree with the corresponding phonon state den- 
sities shown for these metals in Fig.l(f),(g) but a consid- 
erable frequency dependence of the electron-phonon pref- 
actor a^(a;) (indicated by dashed lines) is also predicted. 
For Mo our linear-response calculations are found to rbe 
close to our RMTA calculations and to earlier onesEa. 
The estimated average coupling here is 0.42 which can be 
compared with the value 0.45 deduced from the specific- 
heat measurements (see Table II). The calculated value 



of A for Cu is 0.14. The specific-heat estimate here is 
less reliable possibly because of the smallness of A and 
errors due to the experimental uncertainty in the value 
of 7. Quite likely, however, that there are some errors in 
the DFT value of the density of states connected with 
the many-body effects since copper valence shell 3c?^" is 
close to the strongly correlated 3dP configuration. It 
is known that within DFT the position of d— band is 
higher than experimentally observed. The band which 
crosses the Fermi level is essentially s— band but the ef- 
fect of hybridization with the d— band should lead to low- 
ering the Fermi velocities. The latter effect is stronger 
if the d— band is closer to the Fermi energy. To obtain 
1 + Xs-h 1.1, one has to reduce our calculated N^ep) 
by approximately 20% which is a reasonable estimate 
for the expected influence of the Coulomb correlations. 
Concerning other estimates of A based on the transport 
properties, our calculated Xtr — 0.13 while this value ex- 
tracted from the measured resistivity data is 0.12 (see 
following subsection). Both values are in agreement with 
our superconducting A. 

As the last example, we consider Pd and discuss para- 
magnon effects. The superconductivitj' in Pd is absent 
because of the large spin fluctuationgS. There was also 
a discussion in the literatureEj on the paramagnon con- 
tributions to the mass enhancement in Nb and V. The 
occurrence of paramagnons is connected with the fluctu- 
ations of the electron spins. Paramagnons usually coun- 
teract superconductivity since the latter has its origin in 
the formation of pairs with the opposite spins. To extract 
the paramagnon contribution, we can use out, calculated 
values of A together with the specific-heatEfl estimates 
Xs-h after formula (|l6|). The necessary data are listed 
in Table II. Comparing these results does not leave any 



place for A 



spin 



= Xs 



X in all the materials except in 



Pd which is a typical example for paramagnon effects. 
Here Xs-h = 0.69 and with the use of our calculated 
a'^F{u!) [Fig. 2(h)], the average electron-phonon coupling 
is found to be equal to 0.35. This results in our value of 
Xspi7i — 0.34 for Pd which is close tojts earlier estimate 
0.31 based on the RMTA calculationN. 

After comparing the calculated and experimental spec- 
tral functions, we present the results of our appli- 
cations to solving the Eliashberg gap equation with 
our knowledge of a'^F{uj). Having fixed the Coulomb 
pseudopotential /z*, the superconducting state is now 
completely described by the strong-coupling theory pf 
superconductivitj-ti. According to the Allen-Dynesc£l 
modified McMillancZl formula: 



McM 



exp 



1.04(1 + A) 
A -/!*(! + 0.62A) 



(17) 



the effect of the first reciprocal moment A of a^F(w) on 
Tc is most important. Unfortunately, the estimation of 
the coupling constant from Tc is difficult because of the 
unknown value of . We use standard Matsubara tech- 
nique to solve numerically the Eliashberg equation for Tc 
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and have found /i* which gives the experimental value of 
Tc ■ The cutoff parameters Wcut were taken to be equal 
to ten phonon-boundary frequencies ujmax- To treat the 
Coulomb pseudopotential in terms of the expression (|l^ ) 
when solving the Eliashberg equation we have rescaled 
actually used parameters fi*{ujcut) to /i* — ^*{ujiog) ac- 
cording to 



1 



1 



^J,* fi*{ujcut) 



+ In 



^cut 



(18) 



Table III reports the obtained ^* values. The main con- 
clusion here is that the calculated ^* varies between 0.11 
and 0.17 which is close to the conventional value usually 
taken ^ 0.13. The noticeable exceptions are only Nh and, 
especially, V for which too large /i* have been found. Of 
course, this overestimation occurs over the conventional 
quantity 0.13 while the detailed theoretical data on the 
Coulomb pseudopotential are unknown. We think that 
the obtained quantities are still below the upper limit for 
the allowed fx* values. 

Also listed in Table III are the values of T^'^"^' eval- 
uated after ( p7| ) with our calculated w/og, A and /i*. As 
it can be seen, the exact solution of the Eliashberg equa- 
tion gives Tc (chosen to be the experimental value) which 
slightly deviates from that estimated after the McMillan 
expression. The accuracy of the later is averagely about 
15%. 

To conclude that our spectral functions provide a 
proper description of superconductivity we have found 
the energy-gap parameters Aq. These results are shown 
in the last two rows of Table III. The available tunneling 
data for Al give Aq coinciding with the theoretical value. 
Some overestimation of the coupling constant in compar- 
ing to the experimental one has taken place in Pb, while 
Aq agrees very closely. Let us turn out to the impor- 
tant case of transition metals Nh and Ta. As we have 
discussed already, the main difhculties in the tunneling 
studies in V and Nh are connected with the oxidation of 
the surface layers and the tunneling estimations of the 
coupling constant in Nh vary considerably in the past. 
In contrast to it, the measured superconducting gap in 
Nh has approximately the same salues in all the exper- 
iments and is equal to 1.56 TOeVu. This value perfectly 
agrees with that found by us which is equal to 1.53 meV. 
We have also found a good agreement in the energy gap 
for V which is within 4% of the experimental one. The 
discrepancy in the energy gap for Ta is again 1 — 2%. 
Because of the low transition temperature there is no 
tunneling data for Mo and we only give the theoretical 
value of Ao. We thus see that despite of the discrepancy 
in evaluating the coupling constants, an extremely good 
agreement (1 — 2%) is obtained for the predicted gap 
data. Such a coincidence is readily understood because 
the ratio 2Ao/Tc is slowly varying for different super- 
conductors. (It is 3.52 within the BCS theory.) Fixing 
the Tc to its experimental value makes the value of Aq 
insensitive to the errors in a'^F{uo) and A. 



To summarize, we have found that our results for the 
spectral functions and, in particular, for the coupling con- 
stants are realistic for the correct description of the su- 
perconducting properties. Especially, an excellent agree- 
ment has been found between our calculated A and the 
values extracted from specific-heat measurements. The 
values of A deduced from the available tunnelling exper- 
iments also agree within 10% with our calculations in 
all the materials except Nh and, especially, V. However, 
taking into account the past tendency to correct the tun- 
neling spectra for Nh as well as our calculations for Ta, 
we do not consider these discrepancies as essential. 

d. Transport properties. 

We now report the results of our applications for calcu- 
lating the electron-phonon contribution to the electrical 
and thermal resistivities (conductivities). This field re- 
tains very important and interesting, first, because the 
easily measured transport properties and, especially, the 
electrical resistivities, provide a valuable information on 
the electron-phonon-coupling strength and, second, no 
large-scale investigations of these properties by ab initio 
theoretical calculations appeared so far . 

We calculate the electrical and thermal resistivity us- 
ing the low-order variational approximation (LOVA) and 
our theoretical transport spectral functions found after 
Eq. (|l3|). As follows from Eqs. (|l|) and (|l|) at high 
temperatures: 



nQ,ceiikBT 



^kBN{eF){viy 



(19) 



(20) 



and an important information is contained in the trans- 
port constant Xtr defined by 



duj 



(21) 



It is usually believed that the latter is close to the su- 
perconducting A because the expressions for al^F{ijj) 
and a^F{uj) are quite similar, except for the factor 
[1 — v(k)v(k')/|v(k)p] which preferentially weights the 
backscattering processes. However, there may exist a sig- 
nificant difference between a1^F{uj) aad a^F((jj) for the 
case of strongly nested Fermi surfacea£3 due to the contri- 
bution from the backscattering of electrons between the 
opposite sides of the nested Fermi surface. 

Despite the complexity of the Fermi surfaces in the 
transition metals we have obtained the transport func- 
tions a1j.F{u) quite close to the superconducting a^F(cj) 
which have been shown in Fig.l. The latter is also true 
for the simple metals considered in this work. Unfortu- 
nately, we have not investigated an interesting question 
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about the low-frequency behavior of the afj.F{uj) due to 
a relatively coarse grid of the phonon wave vectors used 
for integrating the Eq. (p^. The values of Xtr calcu- 
lated from our transport functions are listed in Table IV. 
Comparison between Xtr and superconducting constants 
calculated earlier (Table II) gives the difference between 
them within 20% in all the tested materials. This is in 
agreement with previous conclusions that Xtr X for 
transition metalstZI. 

The results of our calculated electrical resistivity p(T) 
and thermal conductivity w~^{T) are presented, respec- 
tively, in Fig. 3 and Fig. 4 (full lines), up to the temper- 
atures 500K. Symbols denote different measured points 
available from Refs. |5^,^. (The residual values of the 
electrical resistivities are subtracted.) For the compari- 
son with the experiment we are limited by the temper- 
atures T < 2Qtr, where Qtr ~ \/ {^'^)tr close to the 
average phonon energy. (We list our calculated values 
of Qtr in Table IV.) This is so because the description 
of the transport properties at high temperatures require 
to take into account the unharmonicity effects and the 
Fermi-surface smearing. At low temperatures (usually 
when T < &tr/5) the calculations will demand the inclu- 
sion of the jAi-sheet and the inelasticity corrections be- 
yond L0VAE3. Also, here a more careful integration over 
the Brillouin zone is necessary to produce a correct limit 
of atrF(uj) when w — > 0. Moreover, we cannot consider 
very low temperatures because of the effects of electron- 
electron scattering, size effects, impurity scattering, etc, 
which may give considerable contributions in addition to 
the electron-phonon scattering. The latter is basically 
responsible for the electrical resistivity of a metal at high 
temperatures in the absence of spin fluctuations. For the 
thermal conductivity, the lattice contribution to the heat 
current also exist and must be taken into account at the 
temperatures at least less than lOOK. From Table IV we 
see that, except lead, the values of Qtr in all other mate- 
rials are well above the low-temperature region and the 
comparison of our results with the measured ones must 
be relevant at the intermediate temperatures. 

To compare the theoretical transport constants with 
the experiment, we fit the measured dataES for p{T) by 
polynomial series 



3-2i 



(22) 



at the temperatures 9(^/2 < T < 2Qtr with n = 2. (The 
accuracy of the fit varies within 3% if n is increased) . The 
empirical values X^r^ were then found using the extracted 
coefficient ci as follows 



where ujp is our calculated bare plasma frequency 
87rN{eF){vl) 



cell 



(23) 



(24) 



The obtained Xtr^ are shown in Table IV where they can 
be compared with our calculated Xtr- Note that similar 
numerical estimates for Xtr^ can also be made by analyz- 
ing thermal conductivities because the measured Lorentz 
number approaches to the Sommerfeld value at the tem- 
peratures T > Qtr- 

For Al [Fig. 3 (a)] we have found a good agreement be- 
tween the theoretical and the experimental resistivities 
at the whole interval of the intermediate temperatures. 
The corresponding values of the transport constants are 
X'if" = 0.37 and A^^^ = 0.39. The reduction of the cou- 
pling constant {Xcaic — 0.44) to the transport one is 
less than 20%. There is also an agreement between the 
theoretical curve and the experimental pointsEHI for the 
thermal conductivity [Fig. 4(a)] above 150i^. The theory, 
however, underestimates w~^{T) at the lower tempera- 
tures. An obvious explanation here is the neglection of 
the lattice contribution to the thermal current. In fact, 
as it can be seen from Fig. 4, such underestimation at the 
low temperatures exist in all other materials considered 
in this work. 

A comparison with the experiment is complicated for 
lead because of its low phonon energies {Qtr ^ 75K) 
and the importance of the unharmonicity effects already 
at the low temperatures. The latter can possibly ex- 
plain our discrepancy in the calculated electrical resistiv- 
ity behavior shown in Fig.3(b). The same disagreement 
exist in our results for the thermal conductivity, Fig. 4(b). 
Here, the discrepancy is minimal at the temperatures 
near 75K and grows fastly as the temperature increases 
mainly because of the linear decay of the measured ther- 
mal conductivity. This obviously contradicts with the 
LOVA behavior of 'w~^{T) and is consistent with our as- 
sumption on the importance of the unharmonicity. Un- 
like in the other considered metals, the computed value 
1.19 of the transport constant here is significantly smaller 
than the electron-phonon A — 1.68. This reduction could 
also point out the importance of the anisotropy in the 
electron-phonon scattering as well as the Fermi-surface 
effects which are not well reproduced by LOVA. In the 
absence of a calculation beyond LOVA it is difficult to 
determine the main source of errors. 

Measured resistivity for Nb starts to saturate at high 
temperatures. In fact, this effect is evident [Fig. 3(d)] at 
the temperatures above 2Qtr ^ AOOK and it does not 
appear at the intermediate interval where the behavior 
of the resistivity only slightly deviates from the LOVA 
prediction. Comparing the calculated Xtr = 1-17 and 
the empirical values A^^^ = 1.11 gives the agreement 
about 5%. Like in Nb, there is a complete agreement 
between the theoretical and the experimental data for V 
[Fig. 3(c)] at the temperatures Qtr/5 < T < 2Qtr (cal- 
culated Qtr ^ 260K ). The theoretical Xtr = 1-15 coin- 
cides with the Xtr^ found empirically. The applicability 
of our method to the description of the transport prop- 
erties for both V and Nb is also supported by comparing 
the thermal-conductivity data. Fig. 4(c), (d). 
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Fig. 3(e) and Fig. 4(e) present the results of our calcu- 
lations in Ta. For w~^{T), an excellent correspondence 
of the theoretical prediction with observed behavior is 
obtained, but p(T) is underestimated in our calcula- 
tion within 10 — 12%. As a consequence, the evaluated 
Xfr^ = 0.93 slightly exceeds the value 0.83 of the theoret- 
ical transport constant. The discrepancy is, in principle, 
within our computaticuial errors. Note however, that the 
experimental behavioiO of p{T) used in evaluating X^^^ 
may depend on the sample purity and should be verified 
by several measurements. The disagreement can also be 
assigned to the influence of high-temperature effects in 
the vicinity of the upper limit {2Qtr ~ MOK) of the in- 
termediate interval . 

The results for Mo are given in Fig. 3(f) and Fig. 4(f). 
To avoid the influence of high-temperature effects we 
have dropped out the measured resistivity points at the 
temperatures above 300-fi' (the calculated Qtr ~ 290K) 
. Fitting for T < 300K gives the result for the empiri- 
cal transport constant A^^^ = 0.40 which is close to our 
prediction, A™''^ = 0.35. The agreement in the thermal- 
conductivity data is satisfactory for the whole interme- 
diate interval. 

The slope of the resistivity in Cu, Fig. 3(g), is also 
obtained quite accurate as in the other materials. The 
value of Xtr is only slightly overestimated in the calcu- 
lation. The discrepancy in the calculated thermal con- 
ductivity. Fig. 4(g), is larger and consists about 20% at 
the room temperature. We cannot explain such disagree- 
ment by the renormalization due to the Coulomb corre- 
lations because, having been proportional to the ratio 
Xtr/i^p, both p{T) and w~^{T) must be insensitive to 
this effect in the first order. The underestimation of 
the theoretical w~^{T) can point out the largeness of 
the lattice contribution to the thermal conductivity at 
the temperatures above lOOK. Unfortunately, there is a 
number of known difficulties to extract the latter values 
from the_experiments£3. From the low-temperature data 
analysisEJ one may conclude that the contribution to the 
thermal resistivity, w^_f^, from the process of a phonon 
decay by emissing electron-hole pairs is very small for Cu 
because of the apparently weak electron-phonon coupling 
(A ~ 0.12 — 0.14). This supports our explanation for the 
obtained discrepancy. 

Finally, the predicted transport properties of Pd are 
presented in Fig. 3(h) and Fig. 4(h). Like Cu, this metal 
has 4cJ^° electronic configuration, but, in contrast to Cu, 
we have found a very good agreement between the calcu- 
lated curve w~^{T) and its measured behavior. We can 
consequently judge that the thermal conductivity carried 
by phonons is pnall in this case which is consistent with 
the conclusionEHl that the contribution w^_f^ is large for 
Pd. We have also found an underestimation of the electri- 
cal resistivity in the calculation. The agreement between 
Xtr and A(^^ (see Table IV) consists about 15% which 
is, in principle, the upper limit of our computational un- 
certainty. Most likely, however, that the additional spin 
fluctuation mechanism of the resistivity is also present in 



this metal. 

In summary, the behavior of p(T) and 'w{T) is consis- 
tent with the results (p^),(p0|) at the intermediate tem- 
peratures and there is no significant discrepancy between 
our calculations and the experimental points. More pre- 
cisely, we have extracted the values of X^^^ using the 
experimental data for p(T) together with our band- 
structure value of Up and found the agreement between 
the experimental and the theoretical transport constants 
to be about 10% (in particular, lower than 5% for 
Al, Nb, Mo and V). In fact, compared with the exper- 
iment is the ratio Xtr/Up. Except possibly Cu and Pd, 
the DFT-based band-structure calculations assumed to 
provide the proper magnitude for the plasma frequency. 
So, we drop the possibility of error cancellations and 
conclude that the theoretical Xtr are in the real agree- 
ment with the experiment. Relatively high error level 
in Pb can be explained by the importance of the unhar- 
monic effects. The large lattice contribution to the ther- 
mal conductivity could affect our comparison for Cu. In 
Pd, the additional mechanism of the resistivity can also 
take place. Nevertheless, taking into account the agree- 
ment in the other calculated properties, we think that our 
description of the electronic transport in the considered 
materials is quite satisfactory. 

IV. CONCLUSION 

We have presented ab initio linear-response calcula- 
tions of the electron-phonon interaction in the transi- 
tion metals Cu, Mo, Nb, Pd, Ta, V and in the sp-metals 
Al, Pb using the local density functional method and the 
LMTO basis set. Our results for the lattice dynami- 
cal, superconducting and transport properties in these 
materials agree well with the experiment. They can be 
summarized as follows: (i) we have obtained tunneling 
spectral functions a'^F{uj) and their first reciprocal mo- 
ments A close to the measured ones; (ii) the correct values 
for the superconducting energy gap have been found us- 
ing our calculated a'^F{uj) and p* corresponding to the 
experimental T^, (iii) the solution of the Eliashberg equa- 
tion for Tc (or for p* if Tc is fixed) is well approximated 
by the conventional McMillan formula; (iv) the mass en- 
hancement observed in the specific-heat measurements 
corresponds very well to our calculations and there is no 
paramagnon contribution in all the metals except Pd; 
(v) we have found the electrical and thermal resistivities 
in agreement with the measured data; (vi) we have also 
found them to be well described by the LOVA expres- 
sions; (vii) the theoretical transport constants agree with 
the values of A^^'' within 10%. To summarize all these 
results, we conclude that our method gives the descrip- 
tion of the electron-phonon coupling with the accuracy 
of order 10%. We also conclude that the effect of renor- 
malization of the energy bands due to electron-electron 
interactions is small in the considered materials. Some 
discrepancies between the theoretical and the tunneling 
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values of A in Nb and V can be assigned to the diffi- 
culty in processing the tunneling data. Nevertheless, it 
seems to us that more experimental and theoretical work 
is necessary to account for the large A and /i* in Nb and, 
especially, in V. 
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TABLE I. Comparison between calculated and experimenta lEj phonon frequencies (THz) at the high-symmetry points X, L 



for the fee metals Al, Pb, Cu, Pd, and at the points H,N for the bcc metals V, Nb,Ta, Mo. Also listed are the theoreti- 
cal-to-experimental volume ratios V/Vq used in the calculations. 



fee |[ bcc 




Al 


Pb 


V 


Nb 


Ta 


Mo 


Cu 


Pd 


Xl \\Hlt 


theory 
exp. 


9.51 
9.69 


1.80 
1.86 


8.03 


6.43 
6.49 


5.13 
5.03 


5.71 
5.52 


7.69 
7.25 


7.17 
6.72 


Xt\\Nl 


theory 
exp. 


5.83 
5.78 


1.06 
0.89 


7.22 


5.52 
5.66 


4.54 
4.35 


7.99 
8.14 


5.36 
5.13 


5.01 
4.64 


Ll\\Nt^ 


theory 
exp. 


9.84 
9.69 


2.18 
2.18 


4.76 


3.94 
3.93 


2.65 
2.63 


5.74 
5.73 


7.77 
7.30 


7.39 
7.02 




theory 
exp. 


4.33 
4.19 


0.92 
0.89 


6.17 


4.80 
5.07 


4.18 
4.35 


4.69 
4.56 


3.64 
3.42 


3.60 
3.34 


V/Vo 




0.955 


1.002 


0.990 


0.972 


0.974 


0.971 


0.985 


0.975 
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TABLE II. Comparison between the calculated electron-phonon coupling constants Xcaic and the values of Xtun deduced 
from the tunnelling experiments. Also listed are the values of Xa^h extracted from the measured specific-heat coefficient 7 with 
the use of our calculated density of states N{eF). 





Al 


Ph 


V 


Nh 


Ta 


Mo 


Cu 


Pd 


Xcalc 


0.44 


1.68 


1.19 


1.26 


0.86 


0.42 


0.14 


0.35 


Xtun 


0.42" 


1.55" 


0.82" 


1.04", 1.22'' 


0.78" 








Xs-h 


0.43 


1.64 


1.00, 1.17 


1.17 


0.83 


0.45 




0.69 



rnJ 
I ' K-'mol 



Ry^cell 



5.49 
1.36"= 



6.87 
3.14^= 



26.14 

9.04^9.82^* 



20.42 
7.66"= 



18.38 
5.84'= 



8.34 
2.10^= 



"Reference \^ 
''Reference 42 
° Reference 54 
''Reference 55 



4.36 
0.69"= 



34.14 
10.0^= 



TABLE III. Calculated values of the Coulomb pseudopotential ^* which provide the experimental values of as the solutions 
of the Eliashberg equation with our knowledge of a^F(ijj). Values T^'''^' were then found with our uJiog, A and /i* in order 
to check the accuracy of the McMillan Tc -expression . Also shown are the computed and the measuredtl superconducting 
energy-gap parameters Aq. 





Al 


Pb 


V 


Nh 


Ta 


Mo 


Cu 


Pd 


* 


0.12 


0.17 


0.30 


0.21 


0.17 


0.14 


0.11 






1.18 
1.22 


7.19 
6.81 


5.40 
6.68 


9.25 
10.5 


4.47 
5.11 


0.92 
0.67 


> 

> 






270 


65 


245 


185 


160 


280 


220 


180 


Ag"'",mey" 

A;;"^mey 


0.18 
0.18 


1.35 
1.33 


0.84 
0.81 


1.53 
1.56 


0.70 
0.71 


0.14 







TABLE IV. Comparison between calculated and empirical values of the transport constant Xtr- The values of X'^^^ were 
deduced from the electrical-resistivity datau with help of our calculated bare plasma frequencies ujp. Also listed are the average 
transport frequencies Qtr = \/ {'-^^)tr 





Al 


Pb 


V 


Nh 


Ta 


Mo 


Cu 


Pd 


\ calc 
Mr 


0.37 


1.19 


1.15 


1.17 


0.83 


0.35 


0.13 


0.43 


\ exp 
\r 


0.39 


1.52 


1.15 


1.11 


0.93 


0.40 


0.12 


0.50 


iJp,eV 


12.29 


14.93 


7.95 


9.47 


9.05 


8.81 


8.75 


7.34 


Qtr,K 


330 


75 


260 


200 


170 


290 


230 


190 
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Fig.l(a)-(h). Calculated phonon-dispersion curves along several symmetry directions for the eight elemental metals 
considered in this work. The lines result frona_the interpolation between the theoretical points (circles) . The results 
of avaiable neutron-difraction measurements^ are shown by triangles. Also plotted are the calculated densities of 
states (DOS). 

Fig.2(a)-(h). Calculated spectral functions a^Flu) of the electron-phonon interaction (full lines) for the eight 
elemental metals considered in this work. The behavior of the electron-pipaon prefactor a^(w) is shown by dashed 
lines. Symbol plots present the results of available tunnelling experimentsOS. 

Fig.3(a)-(h). Calculated temperature dependence of the electrical resistivity, p{T), as a lowest-order variational 
solution of the Boltzmann equation for the eight elemental metals considered in this work. Symbols show different 
experimental data available from Ref. ^ 

Fig. 4(a)-(h). Calculated temperature dependence of the thermal conductivity, w~^{T), as a lowest-order vari- 
ational solution of the Boltzmann equation for the eight elemental metals considered in this work. Symbols show 
different experimental data available from Ref. E^. 
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